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Abstract 



We provide evidence that the KPZ exponents in two-dimensional quantum gravity 
can be interpreted as scaling exponents of correlation functions which are functions 
of the invariant geodesic distance between the fields. 
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1 Introduction 



The calculation of the dressed scaling exponents is a milestone in the theory of 
two-dimensional quantum gravity [0,0. Strictly speaking, however, the derivation 
uses only finite size scaling arguments and involves only matter field correlators 
integrated over all space-time. The result can be formulated as follows: if we are 
given a conformal field theory where a scalar operator has scaling dimension A 
we know that the integrated correlator will scale with the volume V as follows 

d 2 x [ d 2 y (4>(x)4>(y)) ~ V 2 ~ 2A °. (1) 
v Jv 

If the conformal field theory is coupled to quantum gravity the corresponding ex- 
pectation value now includes the average over equivalence classes of metrics with 



/ d 2 xyj g(x) = V and Eq. flU) is replaced by 

(J v d 2 xjgjx) J y d 2 yjg~(y) <t>{x)<j>{y)^ ~ V 2 ~ 2A , (2) 

where the dressed scaling exponent A for a conformal field theory with central charge 
c is related to the scaling exponent A in flat space by 



v 'l -c + 24A /2- ^T=~c 
A ^ 2 ~ v / 25^- VT^c ' ^ 

While these results are beautiful, one key ingredient is missing compared to the 
theory in flat space: the concept of a correlation length which diverges at the critical 
point. In fact, it has occasionally been argued that there exists no such concept in 
quantum gravity since we integrate over all metrics. As we discuss below, however, a 
proper definition of invariant correlation length exists, but even with this definition 
at hand there are situations where the correlation length is not divergent in quantum 
gravity, although it diverges for the same matter system in flat space. For many 
Ising spins coupled to quantum gravity there exist convincing arguments which show 
that the spin-spin correlation length does not diverge at the critical point even if 
the transition is not first order ||. 

The task of this article is to provide evidence that there exists a divergent corre- 
lation length associated to correlation functions of matter fields in two-dimensional 
quantum gravity coupled to matter with central charge in the interval < c < 1. 
First a notation of reparametrization invariant distance is needed in order that one 
can attribute a physical meaning to a divergent correlation length. One trivial sug- 
gestion is the following: for a fixed metric the concept of geodesic distance serves 
well as the definition of invariant distance. We can transport this notation to quan- 
tum gravity provided we perform the quantum average over metrics where the two 
marked points are precisely separated by a geodesic distance R. In this way a 
possible definition of a two-point correlator of <fi{x) could be 

G^R; A) = J V[g]V V e -So-Su J d 2 xd 2 y^ g(x)g(y) 4>(x)cf>(y)5(D g (x, y)-R). (4) 
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In Eq. (|J) D g (x, y) denotes the geodesic distance between x and y with respect to the 
given metric g. The integration T>[g] is intended to be performed over equivalence 
classes of metrics, while T>(p signifies the integration over the matter fields which are 
coupled to quantum gravity via the action Sm(<p)- A is the cosmo logical constant. 

While this definition is straightforward it has not been easy to use in a con- 
tinuum context, i.e. in Liouville theory. An inconclusive attempt was pioneered 
by F. David ||], but the 5- function constraint made it difficult to perform any de- 
tailed calculations. In the context of two-dimensional gravity and non-critical string 
theory, regularized by the use of dynamical triangulations |5|-§, it is possible to 
implement the ^-function constraint in Eq. (^). This was done in ||, using the 
so called transfer matrix formalism of dynamical triangulated surfaces p!0|-p^] (see 



also [13,|14| for continuum transcriptions of this formalism). In |§ the discretized 
version of Eq. (KJ) was calculated in the case = 1 and in the scaling limit the 
following result was obtained: 

Gl (*;A) = A^i^, (5) 
smlr Aii? 

In |l5j this result was generalized to other scaling operators in pure two-dimensional 
quantum gravity and the concept of an operator product expansion was developed, 
the distance of separation between the various operators being geodesic in the same 
way as in Eq. (HI). 

As shown in ||[1(J] an object like Gi(R; A) is a good probe of the fractal structure 



of space-time which is determined by the exponential decay of G\(R; A). We expect 

G 1 {R;A)~ e - cAl/dHR , (6) 

where dn is the Hausdorff dimension of our universes in the following sense: consider 
the ensemble of universes where two marked points are separated a geodesic distance 
R. We define dn by 

(V) R ~ R d » (7) 

provided R is not much larger than l/K l / dH . We conclude from Eq. ([5]) that the 
fractal dimension of space-time in two-dimensional quantum gravity is 4, as first 



proved by different means in |H| 



By a Laplace transformation it is possible to obtain the correlation functions for 
universes of fixed volume V: 



G (j> {R-A) = dV e~ AV G^,(R; V). 

J 



For 0=1 the correlation function Gi(R;V) is related to the average "area" of 
spherical shells Si(R; V) of geodesic radius R by 

where Z(V) denotes the partition function of fixed volume, i.e. 

Z(V) = J V[g\Dif e- Sc - SM 5( J d 2 xjg{x) - V). (10) 
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An alternative and very natural definition of Hausdorff dimension is 

V) ~ R dh ~ l for R < (11) 

By an inverse Laplace transformation of Eq. fl5|) it follows that dh = du = 4. 
However, in the case of matter fields coupled to gravity there exists no proof that 
du = dh- If du = dh it is natural to introduce the finite size scaling variable 

(12) 

From Eqs. ([□]) and (|6]) we expect for Si (R; V) a finite size scaling relation of the 
form 

S^R-V) ~V l - l/dh F(x) , (13) 
where F(x) has the following properties: 

F(x) ~ x^" 1 for x < 1, F(x) -> for x > 1. (14) 

These relations are of course satisfied with = 4 in pure gravity []9,18] and numer- 
ical simulations indicate that it is also true for < c < 1 [I7|, 18]. The finite size 
scaling arguments can also be applied to the more general correlators 

s *<Wtw- (15) 

From Eq. (|3|) we have for a field with scaling exponents A in flat space and A after 



coupling to quantum gravity [17, 18 



dR S^R; V) ~ V 1 -*. (16) 

From this scaling and the definition of dh it is tempting to conjecture the following 
scaling behavior: 

S *( R > ^ ~ fix) = V^-^F.ix), (17) 

where /(0) > and F(x) ~ x d h^-&)-~t f or sma ii Xm This formula generalizes 
Eqs. (|llD-(|14|) to the case where matter is coupled to quantum gravity. The corre- 
sponding formula in flat two-dimensional space reads: 

sf\R;V)-J^f(R/W). (is) 

In this article we report on extensive computer simulations where we tried to 
verify Eq. (|T7j). The finite size aspect of the formula, i.e. the volume dependence, 
has already been checked in ||17|| . However, until now no reliable results exist for the 
correlator itself, i.e. for the R dependence. As explained above this R dependence, 
i.e. the infinite correlation length, is the crucial feature of conformal field theory and 
we provide evidence that it exists also for two-dimensional quantum gravity coupled 
to a conformal field theory if we use the geodesic distance as a measure of length. 
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2 Numerical setup and results 



The numerical simulations are performed as follows: As discretization we use dy- 
namical triangulations. In this formalism surfaces with spherical topology are con- 
structed from equilateral triangles and we allow self energy and tadpole graphs to 
form in the dual (p 3 graph of the triangulation. To each vertex is associated a spin, 
either Ising or three-states Potts spin, depending on the model. A given spin inter- 
acts with the spins on the neighboring vertices. The Monte Carlo updating of the 
triangulations is performed by the so-called flip algorithm and the spins are updated 
by standard cluster algorithms. The flips are organized in "sweeps" which consist of 
approximately Nl accepted flips where Nl is the number of links of the triangulated 
surface. After a sweep we update the spin system. All this is by now standard and 
we refer to Hl9|-|2^| for details about the actions or Monte Carlo procedures. 

The results presented in this paper cover system sizes from 16000 to 128000 
triangles and the number of sweeps is 1.7-5.0 x 10 6 . For the Ising and three-states 
Potts model on dynamically triangulated surfaces the critical temperatures (3 C are 



known |23|,|24|. We use these values in the simulations. 

Geodesic distances on the triangulations are defined either as the shortest link 
distance between two vertices or the shortest path through neighboring triangles. 
While these two distances can vary a lot for individually chosen points on a given 
triangulation, they are proportional when the average over the ensemble of triangu- 
lations is taken. We will report here the results obtained by the use of link distances. 
Our discretized length r will be the number of links. Our discretized volume (or 
area) will be the number N? of triangles used in the computer simulations and our 
discretized scaling variable will be 

where N is the number of vertices on the surface. We denote the discretized distri- 
butions corresponding to S\{R\ V) and S^R; V) by 

ni (r;N) = (E^(Ai-r)), (20a) 
j 

n (r;iV) = (£ a i a j b(D {j - r)) , (20b) 
j 

where the indices % and j label vertices: % is a random fixed vertex for each measure- 
ment and j runs over all vertices of the given configuration. Dij is the link distance 
between the vertices labelled by i and j. 

Our first task is to determine dh- If the scaling hypothesis is correct, both the 
distributions n x (r; N) and n^(r; N) will depend only on the scaling variable x given 
by Eq. (|T§D, except for the overall scaling factors iV 1-1 /^ and jV 1_A-1 M> respectively. 
Hence, we have to determine the best value of dh such that 

n x (r; N) = jV 1_1/dfc F(a:), n^r; N) = N 1 ~ A ~ 1 / dh F^x) (21) 

for all values of r and N. This has already been done for the Ising model and 



the three-states Potts model for the n x (r; N) correlator in [fT^|18ll and for n«^(r; N) 
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in |T7[. Our data is in very good agreement with the results obtained in [T7| ,[ 
namely dh = 4 for both the Ising model and the three-states Potts model coupled 
to gravity Details of these measurements will be published elsewhere [KJ. 



Ising Model on Fixed Lattice c= 3/4 (fixed) G= 4.224(4) B= -4.4(1) 10"(-3) 




r 

W72 



Figure 1: (a) Data for FF lat ( defined in Eq. ( p4[) , for small values of x 
in the case of the Ising model on a fixed triangular lattice with T 2 topology. The 
fit shown is to Eq. ( |25a| ) for N T = 64000 triangles. (6) Same as (a) but the fit is to 
Eq. (25b). 



Let us now turn to the main task, the determination of the short distance behav- 
ior of the spin-spin correlator for the Ising model and the three-states Potts model 
coupled to gravity. The Ising model on a regular lattice has a second order tran- 
sition at a critical temperature /3£ lat and the corresponding conformal field theory 
has central charge c = 1/2. The spin field has scaling dimension A = 1/8 and 
the dressed scaling dimension after coupling to two-dimensional quantum gravity 
is A = 1/3. The tree-states Potts model corresponds to a c = 4/5 conformal field 
theory and the dressed scaling dimension of the spin field in this case is A = 2/5. 
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Ising Model c= 5/3 (fixed) C= 0.945(4) B= 0.0055(2) 



128000 hs^h 
64000 h+-H 
32000 

16000 i ■ 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
X 

Ising Model o= 1 .535(6) C= 0.834(7) 



128000 he-i 

64000 n-^ 
32000 

16000 h« i 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



Figure 2: (a) Same as Fig. p] but for the Ising model coupled to quantum gravity. 
The fit shown is to Eq. (|25aJ) with exponent 4/3 for Nt = 128000 triangles, (b) 
Same as (a) but the fit is now to Eq. (25b). 



The prediction for short distance scaling is then 

F^x) oc z*^ 1 -^- 1 for x < 1, (22) 

i.e. (using the value d h = 4) 

F* sin9 (x) oc x 5/3 for z<l (23a) 
F^ix) oc x 7/5 for (23b) 

In order to calibrate the expected accuracy with which the exponent dh(l— A) — 1 
can be extracted we have performed simulations of an Ising spin system on a flat 
lattice (with periodic boundary conditions) and measured the known function 

F™ at (x) oc x 3/4 for x < 1, (24) 
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Q=3 Potts Model c= 7/5 (fixed) C= 1 .045(1 ) 



128000 hs^h 
64000 h+-H 
32000 
16000 i 



0.05 0.1 0.15 0.2 0.25 0.3 

X 

Q=3 Potts Model c= 1 .396(5) C= 1 .040(8) 



128000 he-i 
64000 
32000 

16000 h« i 



0.2 0.25 0.3 0.35 



Figure 3: (a) Same as Fig. |l|a in the case of the three-states Potts model coupled to 
quantum gravity. The fit shown is to Eq. (|25a|) with exponent 7/5 for Nt = 128000 
triangles, (b) Same as in (a) but the fit is to Eq. (25b). 



where x — rj \[N and r is the distance in lattice units between a spin and its "spher- 
ical shell" . This rather unusual way of measuring correlation functions works quite 
well as one can see from the results shown in Table [I]. We performed simulations 
on regular toroidal triangular lattices of sizes 16000-64000 triangles. In order to 
understand the size of finite size effects and the effect of the function f(x) we show 
the results of the fits to the following functional forms: 

F* la \x) = Cx 3/4 + B, (25a) 
= Cx c , (25b) 
= Cx c + B. (25c) 

The constant B in Eqs. (|25a|) and (|25c|) reflects the fact that the volume element 
is discrete with a smallest unit and the spin correlated with itself in the volume 
element at the shortest distance is 1. The data is shown in Fig. [I] together with the 
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best fit to the largest lattice. It is clear that the determination of the exponent c 
is in excellent agreement with the theory. It is important that we only use a; < 1 
if we fit to Eqs. (p5|). It is pleasant surprise that the continuum formulae are valid 
even for small values of r. For x > 0.06 the function / present in Eq. ( |I8D will play 
an important role. 

For the Ising model and the three-states Potts model coupled to gravity we 
repeat the analysis performed for the regular lattice. We fit the data to Eqs. (25a]), 
with exponents 5/3 and 7/5 respectively, (25b) and ( |25c| ) for x < 0.45. The results 
are shown in Table ^| for the Ising model and Table |3| for the three-states Potts 
model. The corresponding plots including the best fit to the largest lattice are 
shown in Fig. and Fig. |2]& for the Ising model and in Fig. and Fig. |3]& for the 
three-states Potts model. 

Our results are consistent with the exponents conjectured in Eqs. ( |23a| ) and 
(23b). We find that simulating the largest lattices, 128000 triangles, is important 
for obtaining enough data points in the relevant region x < 0.45. The fits to the 
predicted behavior Eq. (|25a| ) are good and the finite size correction B approaches 
zero as the volume increases. Even if the exponent c is allowed to vary, as in the 
fits to Eqs. (25b) and ( |25cj ), it approaches its predicted value convincingly as the 
volume is increased. The small discrepancy is consistent with finite size effects, 
which clearly are more important than the statistical errors quoted in Tables ^| and 
H We find that the difference in the values of c obtained from the fits to Eqs. (25b) 
and ( |25c|) gives a measure of the systematic errors entering from varying the range 
and type of the fits. 



N T 


C B 


c C 


c C B 


16000 
32000 
64000 


4.276(5) -0.0088(2) 
4.248(4) -0.0062(2) 
4.224(4) -0.0044(1) 


0.7806(8) 4.62(2) 
0.769(1) 4.44(2) 
0.765(1) 4.39(2) 


0.754(7) 4.32(7) -0.008(1) 
0.756(4) 4.32(5) -0.005(1) 
0.753(4) 4.25(4) -0.0038(7) 


Table 1: Results for the fits to Eqs. (|25a|), (25b) and (p5c|) for the Ising model on 
a fixed triangular lattice with T 2 topology. 


N T 


C B 


c C 


c C B 


16000 
32000 
64000 
128000 


0.880(4) 0.0173(4) 
0.923(3) 0.0108(2) 
0.940(3) 0.0075(1) 
0.945(4) 0.0054(2) 


1.432(6) 0.770(6) 
1.470(3) 0.793(3) 
1.497(4) 0.805(4) 
1.535(6) 0.834(7) 


1.51(1) 0.799(8) 0.0063(4) 
1.545(5) 0.829(4) 0.0054(1) 
1.565(4) 0.841(4) 0.0043(1) 
1.596(6) 0.870(7) 0.0036(1) 



Table 2: Same as in Table [I] for the Ising model coupled to quantum gravity. 
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N T 


C B 


c C 


c C B 


16000 
32000 
64000 
128000 


0.947(1) 0.0098(1) 
0.993(3) 0.0050(3) 
1.013(2) 0.0024(2) 
1.044(4) 0.0003(3) 


1.291(4) 0.880(5) 
1.348(3) 0.959(5) 
1.371(2) 0.991(4) 
1.396(5) 1.040(8) 


1.331(8) 0.893(6) 0.0056(6) 
1.425(7) 1.007(7) 0.0077(3) 
1.447(5) 1.047(6) 0.0065(2) 
1.470(8) 1.10(1) 0.0055(3) 



Table 3: Same as in Table [TJ for the three-states Potts model coupled to quantum 
gravity. 



3 Discussion 



We have provided substantial evidence for the conjecture put forward in |L8| that the 
the two-point correlator in a unitary conformal field theory translates the two-point 
correlator of the same theory coupled to quantum gravity in the following way: 

Si la \R;V) = ^f fla \R/VV) - S,(R;V) = ^ K f(R/V 1 ^), (26) 

where /(0) > 0, R is the geodesic distance as defined above and dh is the fractal 
dimension of space-time. 

A number of questions still need to be answered before we have a complete 
understanding of the concept of invariant correlation functions in two-dimensional 
quantum gravity. First, the value of dh as a function of the central charge c of 
the conformal field theory should be clarified. Some formal arguments suggest that 
dh = 2m for a (m,m + 1) conformal theory coupled to gravity. This would imply 
that dh —>■ oo for c = 1 — 6/m(m + 1) — > 1. For the three-states Potts model m — 4, 
i.e. the value of dh = 8. This is definitely not seen in the computer simulations, 
where particularly the data for ni(r;N) favor dh = 4. However, one could argue 
that the systems considered so far are much too small for observe dh = 8, since 
( 128000) 1 ^ 8 ~ 4. While this might be true, it contradicts the fact that finite size 
scaling works very well for the three-states Potts model coupled to quantum gravity 
already for considerable smaller systems and we get the correct critical exponents. 
All the scaling arguments of |T7LP| yield dh = 4 with quite small error; the only 



significant discrepancy being that of the height of the maximum of the distribution 
n«^(r; N) for the three-states Potts model. According to Eq. fl2j| ) it should scale 
as i\T d ''( 1-A ) _1 ; but using the theoretical value A = 7/5 we obtain dh = 4.32(2) 
from collapsing the 64A-128A" configurations. Finite size effects, however, play an 
important role since the value of dh extracted this way decreases rapidly with system 
size towards the value 4. We will report more details in the future . 

Next, one should understand the tail of the distribution n^(r; N) for r ^> jV 1//d '\ 
The following simple argument indicates that the distribution should be identical 
to that of pure gravity: for such large values of R the universes we observe are 
essentially one-dimensional, since they have to be long tubes of almost no transverse 
extension. One-dimensional universes cannot (contrary to genuine two-dimensional 
universes) have non-trivial matter interactions unless they are long range, which is 
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not the case here. Consequently there should be no difference between pure gravity 
and gravity including matter fields. This is only a heuristic argument, but if true 
it rules out the existence of a single scaling variable x as in Eq. fllTf ) for the whole 
range of x unless dh = 4, which is the value for pure gravity. The argument might be 
too simple since preliminary results for c = —2 conformal field theories coupled to 
gravity indicate both a common scale variable for all values of x and a dh < 4 |26[ . 
The c = —2 matter theory is on the other hand non-unitary and the heuristic 
arguments might be still be valid for unitary theories. 

Finally, it would be most interesting to be able to prove analytically the results 
reported above and extend the concept of operator product expansion initiated 
in |T5j to the case of unitary conformal field theories coupled to two-dimensional 



quantum gravity. 
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